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A model for the QCD vacuum based on domain-like background gluon fields has previously been developed and 
demonstrated to exhibit confinement of quarks and strong correlation of the local chirality of quark modes and 
duality of the background gluon field. On the boundary of the j-th domain, quark fluctuations satisfy chirality 
violating boundary conditions parametrized by a chiral angle ctj which is treated as a collective variable. The 
dependence of the free energy of an ensemble of N — > oo domains on {aj,j = 1 . . . N} is studied. For one massless 
flavour there are only two degenerate minima {otj = 0, n Vj) due to the contribution of the anomaly to the 
free energy, which suppresses the continuous (7(1) degeneracy. The residual discrete symmetry is spontaneously 
broken with a quark condensate of — (238MeV) 3 . For Nf massless flavours a continuous SU(Nf)L x SU(Nf)a 
degeneracy appears, while the Ua{1) direction remains suppressed. 



1. INTRODUCTION 

Many of the features of quantum chromody- 
namics long expected to emerge in the nonper- 
turbative regime of the theory have been realised 
in a range of phenomenological models as well 
as demonstrated quantitatively in lattice gauge 
theory. Confinement of quarks, spontaneously 
broken flavour chiral symmetry with a nonzero 
quark condensate and resolution of £7a(1) prob- 
lem are believed to be intrinsically connected to 
each other. Models of the QCD vacuum, which 
can give insight into what is going right with 
lattice simulations of QCD, can invariably con- 
nect two of these three phenomena, but rarely 
all three simultaneously. In this paper, we pur- 
sue one step further the exploration of the "do- 
main model" for the vacuum, originally proposed 
in PP, as a scenario for simultaneous appearance 
of all three phenomena: confinement, sponta- 
neous chiral symmetry breaking via the appear- 
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ance of a quark condensate and a continuous 
SU(Nf)z x SU(Nf)n degeneracy of the vacuum 
for Nf massless quarks, but without a ?7a(1) con- 
tinuous degeneracy of ground states that would 
be indicative of an unwanted Goldstone boson. 

The model under discussion is defined by a par- 
tition function describing an ensemble of hyper- 
spherical domains. Each domain is characterised 
by a background covariantly constant self-dual 
or anti-self-dual gluon field of random orienta- 
tion. Summing over all orientations and dualities 
guarantees a Lorentz and CP invariant ensem- 
ble. Confinement of quarks is manifested in the 
model, as demonstrated in the original work pQ. 
On the boundaries of each hypersphere, fermion 
fluctuations satisfy a chirality violating boundary 
condition 

ifj(x)e la ^iP(x) = ip(x) (1) 

which is 2-7T periodic in the chiral angle a. Here ry M 
is a unit radial vector at the boundary. Integrat- 
ing over all such chiral angles guarantees chiral 
invariance of the ensemble. As a consequence of 
Eq. (|TJ the spectrum of eigenvalues A of the Dirac 
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operator in a single domain is asymmetric under 
A — > — A. Such asymmetries have been studied 
in other contexts, for example by [2]. In the case 
of the domain model, the above boundary condi- 
tions are combined with the (anti- )self-dual gluon 
field which leads to a strong correlation between 
the local chirality of quark modes at the centres 
of domains with the duality of the background 
gluon field 3 . In this paper we study how these 
aspects contribute to quark condensate formation 
and the pattern of chiral symmetry breaking in 
the vacuum. 

The vacua of the quantum problem associated 
with an ensemble of domains are the minima of 
the free energy determined from the partition 
function. The problem of the quark contribu- 
tions to the free energy requires calculation of the 
determinant of the Dirac operator in the pres- 
ence of chirality violating boundary conditions. 
For slightly different choice of boundary condi- 
tion (with a — ► — id — it/2) this problem has 
been addressed in 0] where the parity odd part 
of the logarithm of the determinant was identi- 
fied as lndet(« p) ~ 2q$ with q the topological 
charge (not necessarily integer) of the underlying 
gluon field. This is basically the chiral anomaly. 
For the specific gluon field relevant to the domain 
model we have obtained a similar result for the 
parity odd part |5lfi| 



lndet(i p) ~ 2iq(a mod7r), 



(2) 



consistent with 0]. We study the consequences 
of this result in the context of the domain model. 

A nonzero quark condensate is generated in 
the model without there being a continuous chi- 
ral Ua(1) degeneracy of minima of the ensemble 
free energy. There are two factors which inter- 
play to achieve this result. Firstly, when self-dual 
and anti-self-dual configurations are summed, the 
anomaly Eq. J2J leads to a contribution to the 
free energy of the form — lncos(2<7arctan(tana)) 
which vanishes when a = 0, 7T. The vacuum is 
degenerate with respect to -Z2 chiral transforma- 
tions. The second factor is the spectral asym- 
metry itself which in the presence of an infinitesi- 
mally small quark mass leads to an addition to the 
free energy linear in the mass, which removes de- 
generacy between two discrete minima, and thus 



generates the nonzero quark condensate. This 
gives a model with the chiral Zi discrete sub- 
group of C/a(1) being spontaneously broken, and 
not the continuous Ua(1) itself. In the absence 
of the mass term the ensemble average of ipip cor- 
rectly vanishes. 

Moreover, the form of Eq.0 means that the 
free energy does not depend on flavour nonsinglet 
chiral angles when more than one massless quark 
flavours are introduced. This allows for the cor- 
rect degeneracy of vacua with respect to continu- 
ous SU(Nf)z, x SU(Nf)n chiral transformations. 
This vacuum structure implies the existence of 
Goldstone bosons in the flavour nonsinglet pseu- 
doscalar channel but not in the singlet channel, 
which a study of the structure of pseudoscalar 
correlation functions unveils. 

2. THE DOMAIN MODEL 

For motivation and a detailed description of the 
model we refer the reader to fj]. The essential 
definition of the model is given in terms of the 
following partition function for N — ► 00 domains 
of radius R 



Z = N lim TT / da, f 
V,Ar^oo- L - L / /r. 



VQ l S[D(B^)Q^]A FP [B^\Q 



Wl 



xe 



(3) 



where the functional spaces of integration Tq and 
T\ are specified by the boundary conditions (x — 

zC) 2 = R 2 



ftiQ®(x)=0, 



(4) 
(5) 
(6) 



Here hi = n°;t a with the generators t a of SU C (3) 
in the adjoint representation and the a,i are chi- 
ral angles associated with for the boundary con- 
dition Eq.© with different values randomly as- 
signed to domains. We shall discuss this con- 
straint in detail in later sections. The thermo- 
dynamic limit assumes V, N — *■ 00 but with the 
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density v^ 1 — N/V taken fixed and finite. The 
partition function is formulated in a background 
field gauge with respect to the domain mean field, 
which is approximated inside and on the bound- 
aries of the domains by a covariantly constant 
(anti-)self-dual gluon field with the field-strength 
tensor of the form 

N 
3=1 

B$B$ = B 2 S up . (7) 

Here are the positions of the centres of domains 
in Euclidean space. 

The measure of integration over parameters 
characterising domains is 

= 48^X^1 dai 

s 

//■7T /■ Z7T 

difii / d-di sin #j / d£i 
Jo Jo 

o 

^ , (p (2Z + 1)7T 

x 2^ 5(6 — ) 

(=0,1,2 

x / dui S(uji - nk) . . . , (8) 
70 fe=0,l 

where are the spherical angles of the chro- 

momagnetic field, u>i is the angle between chro- 
moelectric and chromomagnetic fields and & is an 
angle parametrising the colour orientation. 

This partition function describes a statistical 
system of these domain-like structures, of density 
v^ 1 where the volume of a domain is v = 7r 2 _R 4 /2, 
and where each domain is characterised by a set of 
internal parameters and whose internal dynamics 
are represented by fluctuation fields. It respects 
all the symmetries of the QCD Lagrangian, since 
the statistical ensemble is invariant under space- 
time and colour gauge symmetries. For the same 
reason, if the quarks are massless then the chiral 
invariance is respected. 

Thus the model involves only two free param- 
eters: the mean field strength B and the mean 
domain radius R. Within this framework the 
gluon condensate to lowest order in fluctuations 



is AB 2 and the topological charge per domain is 
q = B 2 R A /16. More significantly, an area law is 
obtained for static quarks. Computation of the 
Wilson loop for a circular contour of a large ra- 
dius L 3> R gives a string tension a = Bf(irBR 2 ) 
where / is given for colour SU(2) and SU (3) in 
and has a purely geometrical meaning in terms of 
overlaps of hyperspheres. Estimations of the val- 
ues of these quantities are known from lattice cal- 
culation or phenomenological approaches and can 
be used to fit B and R. As described in ^ these 
parameters are fixed to be y/~B = 947MeV, R = 
(760MeV) _1 = 0.26fm with the average absolute 
value of topological charge per domain turning 
out to be q w 0.15 and the density of domains 
u _1 = 42fm~ 4 . The topological susceptibility 
then turns out to be x ~ (197MeV) 4 , comparable 
to the Witten-Veneziano value [7J- This fixing of 
the parameters of the model remains unchanged 
in this investigation of the quark sector. The 
quark condensate at the origin of a domain where 
angular dependence drops out was estimated in 
paper [T] with a result -(228MeV) 3 . 

3. DIRAC OPERATOR AND SPEC- 
TRUM 

The eigenvalue problem 

W>(x) = XtP(x), 
if,{x)e ta ^^{x) = iP(x),x 2 =R 2 (9) 

was studied in Dirac matrices are in anti- 
hermitean representation. For a assumed to be 
real a bi-orthogonal basis has to be constructed. 
Solutions can be labelled via the Casimirs and 
eigenvalues 

K 2 = K2-£(| + l), fe = 0,l,...,co 

K l,2 m l,2, 

mi, 2 = -fc/2,-fc/2 + l,...,fc/2- l,k/2, 

corresponding to the angular momentum opera- 
tors 

K 1)2 = i(L±M) 

with L the usual three-dimensional angular mo- 
mentum operator and M the Euclidean version 
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of the boost operator. The solutions for the self- 
dual background field are then 



(10) 



where x an d f must both have negative chirality 
in the self-dual field and k related to the polari- 
sation of the field defined via the projector 



O k = N + E K +N.E_ K 
with 

N± = ~(l±n/|n|), E d 



(11) 



-(1±SB/B) 



being respectively separate projectors for colour 
and spin polarizations. Significantly, the nega- 
tive chirality for \ an d V is the only choice for 
which the boundary condition Eq.JSJ) can be im- 
plemented for the self-dual background. The ex- 
plicit form of the spinors x an d f can be found 
in pi], where it is demonstrated that the eigen- 
spinor Eq. (|10|l has definite chirality at the centre 
of domain correlated with the duality of the gluon 
field. The boundary condition reduces to 



X 



-e Tl >, X 



ipe 



R 2 



(12) 



where upper (lower) signs correspond to <p and 
X with chirality =pl, which, using the solutions, 
amounts to equations for the two possible polar- 



isations, for A fc + : 



e' la M (k 

x [M (k + 2 
k + 2-A 2 



2- A 2 ,k + 2,z ) 
- A 2 ,k + 2,z ) 



iA 



k + 2 



M(k + 3- A 2 ,/s + 3,z ) 



= 0, 
and for A 



k 

Mi 



-A 2 ,k + 2 7 z ) 



+ , v V /(l-A 2 ,fc 



k + 2 



3, z ) = 



(13) 



(14) 



where z = BR 2 /2 and A = X/V2B. For the 
present work Eas. (|13ll4fl are the starting point, 
from which we see by inspection that a dis- 
crete spectrum of complex eigenvalues emerges 
for which there is no symmetry of the form A — > 



— A. For given chirality and polarisation and an- 
gular momentum k, an infinite set of discrete A 
are obtained labelled by a "principal quantum 
number" n. 

4. QUARK DETERMINANT AND 
FREE ENERGY FOR A SINGLE DO- 
MAIN 

We consider the one-loop contribution of the 
quarks to the free energy density F(B,R\a) of 
a single (anti-)self-dual domain of volume v = 

Tr 2 Ry2 



exp{~vF(B,R\a)} = det Q 

= n 



/t,fc,n,mi 

= exp{-C'( S )} s=0 . (15) 

The normalization is chosen such that 
lims-fO F(B, R\a) — 0. The free energy is then 
F = u _1 C'(0). In [516) . we discuss in more detail 
the technical aspects of the computation of this 
quantity in zcta function regularisation, start- 
ing with a purely imaginary choice for a which 
provides for a real spectrum, and analytically 
continuing the final result to real values of a. It 
is known that £(s) has two contributions: from 
Qj2 (s), which is insensitive to the spectral asym- 
metry, and from the spectral function ry(s) which 
is odd under A — > —A. These two functions are 
defined by 



'(*) = 



k.n.K 



(fc + i) 



2 s 



[AL(0)P 

V (s) = /i'£(fc + l) 



k,n,K, 

sgn(Ag n (0)) 



Bgn(Ag n (B)) 

\Kn(B)\ S 



(16) 



(17) 



In terms of these, the free energy density for a 
given parameter a is: 



F = V- 1 Q<^(0)±i^: 



<0) + l - V (0) 



(18) 
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with /i - arbitrary scale. The final result for F is 
complex with the imaginary part of the form 



Q?F = ±2gArctan(tan(a)) 



(19) 



where q is the absolute value of topological charge 
in a domain. This charge is not integer here in 
general but the anomalous term is im periodic 
in a. This is the Abelian anomaly as observed 
within the context of bag-like boundary condi- 
tions by 0] . Its appearance here is in the spirit of 
the derivation by Fujikawa [Bj, where the phase 
appears as an extra contribution under a chiral 
transformation on the fermionic measure of inte- 
gration. However, in our calculations |5I6| . an a 
dependent real part also appears. To explain its 
possible fate some technical details of the calcu- 
lation are necessary. The sum over modes n is 
converted into a contour integral via the relation 

sis £ £!■■"«> (20 > 

where the zeroes of / define the eigenvalues. Thus 
the left hand side of the constraint equations 
Eas. l|13ll4|l appear in these expressions. One 
then performs a Debye-like expansion in l/k of 
the Kummer functions. Exchanging this summa- 
tion with the sum over k one can read off each 
order in l/k via an ordinary Riemann function 
whose analytic behaviour is known. However in 
principle the order of summations does not com- 
mute - there can be extra terms, which are diffi- 
cult for computation even in the simplest exam- 
ples 9 . Taking into account such contributions 
is expected to remove the a dependence from 
the real part, leaving only the anomaly contri- 
bution from the imaginary part, consistent with 
0] . Thus for the purpose of the following investi- 
gation, we will assume that the anomaly Eq. (|19J) 
provides the entire result for the a-dependent 
part of the free energy of massless fermions in 
a domain. 

The chiral condensate is computed via the pres- 
ence in the free energy of a term linear in an in- 
finitcsimally small mass, namely 2 



F = F m = + i — ry(l). 



(21) 



Using the analogue of the representation Eq. (|20|l 
the summed contributions of both polarisations in 



the self-dual domain to the asymmetry function 
can be written, 



r](s) = i[iRe l 



, cos(tts/2) v^ /; ,i- s 

7r(l - s) i 

y ' k=l 

x [1 + M(l,k + 2,z) 



M(l,fc + 3,-z) 



fc + 2 



(22) 



We next evaluate the asymptotic behaviour in k. 
A singular term as s — > 1 can be extracted, which 
turns out to be field (that is, B) independent and 
is canceled by the normalization. The final ex- 
pression for the term linear in mass in the free 
energy density for a self-dual domain is 

Spied) = _ e « TO 

where we have used the suggestive notation 
1 ^ k 



ir 2 R 3 ^ k 
x [M(l,fc + 2,z) 



1 



fc + 2 



M(l,fc + 3,-z) - 1 



(23) 



coming from 77(1) with the sum over z correpond- 
ing to a color trace. The a dependent part of 
the free energy of a self-dual domain for mas- 
sive quarks is complex with the following real and 
imaginary parts: 



F 

5RF 



= SJF + i%F 



-to cos a -< iptf> y 



= 2-Arctan(tan(a)) 

v 

— to sin (x -< ipif] >- . 



(24) 



(25) 



The free energy of an anti-self-dual domain is ob- 
tained via complex conjugation. 

5. GROUND STATE AND CONDEN- 
SATE 

Under the assumption that only the anoma- 
lous term depends on the chiral angle, the part of 
the free energy density T relevant for the present 
consideration of an ensemble of N — > 00 domains 
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with both self-dual and anti-self-dual configura- 
tions takes the form 



-vNF 



N ,2; 



= da, 

= A/"exp ^iVmaxln(cos(w3i 7 '(Q!)))^ . 

The maxima (minima of the free energy density) 
are achieved at ai = . . . = a at = ^n. 

In the absence of a quark mass, only the 
anomaly contribution in the imaginary part, 
S-F, of the free energy of a single domain ap- 
pears under the logarithm of the cosine and 
defines the minima of the free energy density, 
ln(cos(w3F(a))) = 0. Thus for massless quarks 
there is no continuous Ua{^) symmetry in the 
ground states rather a discrete Z 2 chiral sym- 
metry The anomaly plays a peculiar role here: 
selecting out those chiral angles which minimise 
the free energy so that the full ?7a(1) group is 
no longer reflected in the vacuum degeneracy. It 
should be stressed here that this residual discrete 
degeneracy is sufficient to ensure zero value for 
the quark condensate in the absence of mass term 
or some other external chirality violating sources. 

In a different and more general context the 
idea of dynamical relaxation of the effective 6- 
parameter based on an energy minimization cri- 
terion was discussed in great detail in |1U| . 

Now switching on the quark mass, we see this 
discrete symmetry spontaneously broken, and one 
of the two vacua selected in the infinite volume 
limit according to the sign of the mass. In this 
case (for these conventions of boundary condition 
and mass term), it is the minimum at a = 0. The 
quark condensate can be now extracted from the 
free energy via 



(ip(x)ip{x)) 



lim lim (vN)~ 

m— >0 N— *oo 



dm 



then m — > gives a nonzero condensate 

(ip(x)ip(x)} = — < ipip y . 

According to Eq. H23|) the condensate is equal to 

(i>{x)ip(x)) = -(237.8 MeV) 3 

for the values of field strength B and domain ra- 
dius R fixed earlier by consideration of the pure 
gluonic characteristics of the vacuum - string ten- 
sion, topological succeptibility and gluon conden- 
sate. A nonzero condensate is generated without 
a continuous degeneracy of the ground states of 
the system. 

6. MULTIFLAVOUR CASE 

The question remains though whether any con- 
tinuous directions in the vacua to be expected 
when the full flavour chiral symmetry is brought 
into play. For this we must generalise the anal- 
ysis. We consider Nf massless quark flavours. 
Firstly, we observe that the fermion boundary 
condition in Eq.J^J), explicitly breaks all chiral 
symmetries, flavour singlet and non-singlet (see 
also 0]). Thus the procedure we have used here 
of integrating over all a does not suffice to restore 
the full chiral symmetry of the massless QCD ac- 
tion. Rather, the boundary condition must be 
generalised to include flavour non-singlet angles, 



a - 



(3 a T° 



(26) 



with T a the Nj - 1 generators of SU(N f ). 



Nj angles 



Taking the thermodynamic limit N — > oo and 



Then integration over 1 + Nj — 1 
a, (3 a S [0, 2n] must be performed for a fully chi 
ral symmetric ensemble. The spectrum of the 
Dirac problem now proceeds quite analogously, 
except that the boundary condition mixes flavour 
components, thus an additional projection into 
flavour sectors is required in order to extract the 
eigenvalue equation analogous to Ea. l(T2^l . 

For Nf = 2 the projection is easy to write 
down, (1 ± f3 ■ tj)/2, which will have the effect 
that a in the dependence of eigenvalues and spec- 
tral functions is replaced by a± \(3\/2 for the two 
Nf = 2 isospin projections. Thus for a given do- 
main, the imaginary part of the free energy will 
be the sum of two terms involving the imaginary 
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part determined for a single flavour but evaluated 
at the two shifted angles: 

3F(a + \0\/2) + - 0\/2) = $SF(a) (27) 

since ^SF(a) — 2q arctan(tan(a)) = 2q(a + rm). 
This expresses the known fact that the anomaly 
only appears in the flavour singlet sector or is 
Abelian. Thus for an ensemble of domains the 
free energy is identical to that for one massless 
flavour, namely it depends only on the Abelian 
angle a. Thus for Nt = 2, the Ua(1) direction 
remains fixed by energy minimisation while the 
SU(2)l x SU(2)r directions represent degenera- 
cies in the space of ground states in the thermo- 
dynamic limit. 

The generalisation to SU (Nf) l X SU (Nj) R can 
be summarised as follows: the functions ±|/3| in 
Eq. i|27[) become functions B i ((3 a ) which corre- 
spond to the number of diagonal generators of the 
flavour group. The cancellation in Ea. (|27l) gener- 
alises reflecting the tracelessness of the SU(Nf) 

generators, Y.i B KP a ) = °- Thus for N f = 3 
one expects eight not nine continuous directions 
in the space of vacua. 

The consequence of this realisation of chiral 
symmetries in the meson spectrum can be seen in 
the general structure of correlators of the flavour 
nonsinglet Jp(x) and singlet Jp(x) pseudoscalar 
quark currents as they appear in the domain 
model, 

(J£(x)4(y)) = gggggj 
(Jp(x)Jp(y)) = ((Jp(x)Jp(y))) 

- «Jp(x)»«Jp(2/)». 

Here double brackets denote integration over 
quantum fluctuation fields and the overline means 
integration over all configurations in the domain 
ensemble. The second line in the right hand side 
of the flavour singlet correlator reflects the inho- 
mogeneity of the background field for a particular 
domain configuration, translation invariance be- 
ing restored only for the ensemble average. The 
second term is entirely determined by the cor- 
relation function of the background gluon field B 
in the ensemble pQ . The analogous "disconnected 
term" in the flavour nonsinglet correlator is equal 
to zero due to the trace over flavour indices. It 



should be added that the pseudoscalar conden- 
sate, (ip"f5ip}, naturally vanishes since parity is 
not broken in the ensemble of domains. Thus 
massless modes can be expected in the nonsin- 
glet channel, but not in the flavour singlet due to 
the additional term in the correlator. This gen- 
eral structure of correlators is exactly the same 
as in the instanton liquid model 1 1 1 j but mani- 
fests the mechanism for eta-prime mass genera- 
tion proposed by Witten in |7| and appreciated 
in chiral perturbation theory by |12|. 
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